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Abstract 

This paper studies the combinatoric structure of the set of all representations, up to equivalence, 
of a finite-dimensional semisimple Lie algebra. This has intrinsic interest as a previously unsolved 
problem in representation theory, and also has applications to the understanding of quantum 
decoherence. We prove that for Hilbert spaces of sufficiently high dimension, decoherence-free 
subspaces exist for almost all representations of the error algebra. For decoherence-free subsystems, 
we plot the function fd{n) which is the fraction of all d-dimensional quantum systems which 
preserve n bits of information through DF subsystems, and note that this function fits an inverse 
beta distribution. The mathematical tools which arise include techniques from classical number 
theory. 

PACS numbers: 03.65.Yz, 03.67.Pp, 03.65.Fd, 02.20.Qs 
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I. INTRODUCTION 



quantum error correction codes 
subsystems jl^ , bang-bang decoupling 



A variety of schemes for protectjii^ quantum information have been developed, including 

decoherence free subspaces jsj O, Q, O, S| , noiseless 
I, and topological quantum computation [l^. The 
first four of these techniques are closely related to each other and can be described in a 
simple unified framework based on representations of the algebra of errors jiol. Iisl ll^ . More 
recently, it was shown [15] that topological quantum computation is also related to the 
error-algebra framework. 

Decoherence-free subspaces and subsystems have already been realized in many impor- 
tant laboratory experiments, of which we mention a few. Consider a system exhibiting 

« I 

electromagnetically induced transparency jilg], such as a tenuous vapor of 3-level atoms in 
ihe A configuration (e.g. Strontium), contained in an optical resonator, as described in 



17[. Decoherence in this system arises from spontaneous emission, which corresponds to a 
transition from the excited state to any of the lower lying states. Since spontaneous emis- 
sion occurs only from the excited state, those states in the orthogonal complement are a 
decoherence-free subspace. 

Zanardi and Rasetti Q] showed that decoherence-free subspaces exist in the phenomeno- 
logically important spin-boson model of spins coupled to a bosonic quantum field, under 
the assumption that the coupling constants are the same for several of the spins. DFSs 
also arise as a consequence of a collective coupling to the environment in the Jaynes- 
Cummings Hamiltonian for A^ identical two-level atoms coupled to a single mode radiation 

field BQQ. 

In this short note, we consider all possible representations of the error algebra 5U3, and 
determine how many of these admit decoherence-free subspaces and noiseless subsystems of 
various sizes. We prove that for quantum mechanical systems of sufficiently high dimen- 
sion, decoherence-free subspaces exist for almost all representations of the error algebra. 
Therefore, our results provide further arguments in favor of the error algebra framework. 

Consider a quantum mechanical system with a finite-dimensional Hilbert space Ti, and 
let D = dim(H). This class of systems includes all spin systems, and the state space of any 
quantum computer. Further suppose that our system is an open quantum system, which 
means that it interacts with another system, called the "environment," and the latter has 
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a Hilbert space 8. Unitary evolution for the combined system Ti. 8 implies non-unitary 
evolution for the reduced density matrix p : H ^ H obtained by tracing out the degrees of 
freedom associated to 8. 

The interaction Hamiltonian H : 1-L®8 ^ 1-L®8 which describes the system-environment 
coupling necessarily takes the form 

k 

H = Y,Fa®Ba, (1) 

a=l 

where {Fa : a = 1 ... A;} is a set of Hermitian operators on Ti. 

If the commutators [-F^, F^] are expressible as linear combinations of the F^, 

[Fa,F,]=lY,fab'F,, (2) 

c 

then according to a classic result on decoherence Q], a subspace V dTi is decoherence-free 
if and only if V decomposes as a direct sum of singlets under the action of the Fa operators, 
where a singlet is defined to be a 1-dimensional invariant subspace of V . 

Assuming 0, we may define an abstract Lie algebra q with generators {xa}, where 
a = 1 . . .k and relations [xa, Xb\ = ifab^Xc- There is then a unique linear map 

<Pf:Q^ End(7^) (3) 

with 0_F(xa) = Fa- The map (^p is a unitary representation of the Lie algebra, since the Fa 
are Hermitian by construction. The reason we make a point to distinguish the abstract Lie 
algebra g from its generators in a particular representation is that we wish to consider the 
effect of changing the representation without changing the symmetry algebra 0. 

For some dimension-dependent number p between and 100, decoherence-free subspaces 
exist in exactly p percent of all Z)-dimensional open quantum systems with errors acting in 
all possible ways. In order to allow concrete calculations, we restrict attention to the error 
algebra g = SU3; the calculations for other Lie algebras differ only in technical aspects. 

If p < 5%, one might be tempted to question the utility of the error-algebra framework, 
since representations which satisfy the necessary criteria would seem to be sparse in the set 
of all representations. On the other hand, if p is close to 100%, then reliable methods for 
protecting quantum information from su(3)-generated errors are plentiful. Therefore, the 
precise value of p is of fundamental importance. In this paper, we show that a reasonable 
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answer for D on the order of 10^ ~ 10^ is j9 ~ 70-90 percent, and p increases with D. Also, 
p seems to asymptotically approach 100%, so in higher dimensional systems there are many 
more representations which protect quantum information (provided that the error algebra 
does not also increase in dimension). 

Our results imply that in most dimensions, decoherence-free subspaces exist for almost 
all representations of the error algebra. However, there are certain dimensions to watch out 
for: in D = 6, only three of eight total representations contain a singlet, or 37.5 percent. 
Therefore, it's especially difficult to preserve quantum information in a 6-dimensional Hilbert 
space, when the errors are generated by SU3. Interestingly, this is related to the fact that 3 
divides 6. 

While our results are mathematically rigorous, any connection of these results to labo- 
ratory studies of decoherence must be classified as highly speculative, for reasons which we 
now explain. A finite-dimensional open quantum system is specified by a Hilbert space Ti, 
system Hamiltonian Hs G End(H), error algebra g, and a representation 0:0^ End(7-^). 
This paper merely determines (for fixed 7i, iJ^, g) the number of representations (f) that sup- 
port protection of quantum coherence, and the statistical distribution of the sizes of noiseless 
subsystems over the space of possible representations of the error algebra. We do not dis- 
cuss the probability that any particular physical system evolves according to the dynamics 
specified by such a representation. 

As a thought experiment, consider an initially isolated D-level quantum system, which is 
then brought into contact with an environment in a controlled way, so that (we suppose) the 
error algebra is known to be SU3. One would like to know the true likelihood of encountering 
a decoherence-free subspace in this situation. The numbers reported in Fig. |21 would answer 
this question if, by some unknown mechanism. Nature selects one of the available 5U3- 
representations at random, with a uniform weighting. As no mechanism which operates this 
way is currently known, it seems impossible at present to connect Fig. |2l with experiment. 

On the other hand, our results show where to look in order to find examples with noiseless 
subsystems, and also highlight situations in which DF subspaces or noiseless subsystems are 
rare. As such, these results should help theorists to construct models in which quantum 
coherence is protected. 
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II. AN EXACT FORMULA FOR THE NUMBER OF SU(3) REPRESENTATIONS 
IN A FIXED DIMENSION 



A. The Exact Formula 

Let ^(n) equal the number of irreducible representations of the Lie algebra su(3) in 
dimension n. We compute an explicit formula for this function in Appendix 1X1 

We now fix a dimension D and determine certain facts about the structure of all possible 
representations (including reducible ones) in that dimension. In particular, we determine the 
fraction which contain a singlet. Any representation in dimension D may be decomposed into 
irreducibles, and the dimensions Pi of those irreducible components determine a partition 
[20I 1 of the integer D. We say that the representation has shape corresponding to a given 
partition if the dimensions of its irreducible components are precisely the integers appearing 
in that partition. 

For a partition p = (pi, . . . ,p„) of D, a naive guess for the number of representations 
with shape p would be nr=i^fe)- '^^^^ guess is correct if and only if p does not contain 
repetitions. For example, if pi = p2 and C,{pi) > 1 then we are over-counting. In order to 
properly count these cases, we use a standard combinatorial function, which we now explain. 

If Pi is repeated rii times in a given partition, we think of the k = C,{pi) different pi 
dimensional irreps as "letters" in an alphabet. The set of distinct n^-fold direct sums of 
these irreps is in 1-1 correspondence with the set of length rii words in k letters, where order 
of letters in a word is not important. 

Lemma 1. Let S{n,k) equal the number of length n strings from an alphabet of k letters, 
with order not important. Then 



The recursion relation S{n, k) = Yll=o '^'(^ — i,k — 1) follows immediately, and gives an 
efficient way of calculating the values of 5*. We now return to our objective of counting the 
total number of su(3) modules in dimension D. 

Let V{D) denote the set of all partitions of D. Given a set H of non-negative integers, let 
V{H, n) denote the set of partitions n = Yli.Pi with pi ^ H for all i, so that V{H, n) C V{n) 
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for all if C N. Let R3 be the set of possible dimensions of an sus-module, 



R3 = {deN: ^ 0} 



{1,3,6,8,10,15,21,...}. 



(4) 



Therefore computation of reduces to computation of ^, which is done in the appendix. 

Theorem 1. Let Mod{su3, D) denote the total number o/su(3) modules in dimension D. 
Assume that each partition p of D is to be expressed in the form D = Yl'i=i "^iPi where the 
Pi are all distinct. Then 



When ^ 7^ 0, its most likely value is = 2, and then the binomial coefficient simplifies to 
(n + 1 choose n) = n + 1. In that case, the partitions that contribute the most are those 
which maximize the product HiLil^* + -'-)' which is the same as maximizing Y\.^i- Therefore, 
the largest terms in the sum are those that do not contain singlets; however, the terms 
which do contain singlets are more numerous. The competition between these two types of 
terms determines the fraction of representations which contain a singlet, which we analyze 
in the next section. 

B. An algorithm for efficient computation 

The sum in (jS} is over V{Rs, D), the set of partitions of D with parts in R3, a specific 
subset of the positive integers. For efficient computation of (0), it is essential to have an 
algorithm which lists only the partitions we are interested in, without having to first list 
all partitions and then filter out those which do not meet our criteria. As the number of 
partitions grows exponentially according to the Hardy- Ramanuj an formula, any reduction 
in the number of terms is crucial to make computation of Mod(su3, n) even possible. 

Fix a set of non-negative integers. 



with Hi < Uj if i < j and ni > 1. We present a simple algorithm for explicitly computing 
the set of partitions of n whose parts lie in H. 





H = {l,ni,n2,...} gN 
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A partition n = YllLi is said to be in reverse lexicographic order if ki > k2 >■■■> km- 
We will refer to a partition which is in reverse lexicographic order as an ordered partition 
for brevity. 

Let P{n, k) compute all ordered partitions of n that begin with a number between 1 and 
k. Let 



denote the subset consisting of those ordered partitions whose elements come from H. Sim- 
ilarly define Pii{n). 

If K = {fci, . . . , kp} is an ordered partition of an integer fc, and £ > is a positive integer, 
let £ V K denote the ordered partition oi k + I given by prepending £ to k, i.e. 



We extend this notation to sets of partitions in the obvious way, so that d V P{n) denotes 
V K : K G P{n)}. 

Theorem 2. The set of partitions whose parts come from H is given by 



Theorem 121 gives a fast recursive algorithm for enumerating only the partitions of n with 
parts in H. Even on a fast computer, the naive algorithm of enumerating all partitions and 
subsequently choosing those with the right properties will lead to space problems due to the 
exponential growth of ©• 

C. Asymptotics for large n 

Formula © gives a way of evaluating Mod(5U3,n) numerically, although this is not effi- 
cient (or even possible at all for large n) as it involves a sum over partitions. If p{n) denotes 
the number of all partitions for n, the Hardy- Ramanuj an asymptotic formula gives 



P(n) = [jP{n,k) 



denote all ordered partitions. Let 



Pnin, k) = {{ki, ...,km)e Pin, k) \ h e H Vz} 



^ V K 



{£, fci, . . . , kp]. 




p{n) 



1 



exp TT 



(6) 
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Good references on this include 20L |21|. 



Conjecture 1. Mod(su3,n) is given, asymptotically for large n, by a formula similar to the 
Hardy-Ramanujan formula (jH)) . Specifically, 

Mod(su3,n) ~ —exp(bn'^) 
n 

for some positive, real constants a, b, c. 

The inspiration for this conjecture was Meinardus' theorem (see [3] and references 
therein), which states roughly that a large class of partition functions have exponential 
behavior which generalizes that of the Hardy-Ramanujan function. 

Nonlinear regression over the interval [1, 110] determined that 

(a, b, c) ~ (0.0771591, 2.70605, 0.459802) (7) 

which is surprisingly close to the values 0.14, 2.57, and 0.5 given by ©. Fig. ^ shows the 
curve (a/n) exp(6n^) with the approximate values ((Tj), together with points representing the 
exact values of Mod(su(3),n) computed using Theorems [T] and |2l for n = 1 (mod 3). 

III. STATISTICS OF THE SET OF ALL REPRESENTATIONS 
A. Singlets 

As an example to illustrate formula consider D = 6. There are 11 partitions of 6, 
only 4 of which correspond to direct sum decompositions of representations of SU3: 

p = {6} 2 representations 

p = {3, 3} 3 representations 

p = {3, 1, 1, 1} 2 representations 

p = {1, . . . ,1} 1 representation 

Therefore in D = 6, three of eight total representations contain a singlet, or 37.5 % . By 
contrast, for D = 5, because there are no representations of su(3) in dimensions 2, 4 and 5, 
only two partitions contribute. These are 5 = 3 + 1 + 1 and 5 = 1 + ■ ■ ■ + 1 giving 2 and 1 
representations respectively. Thus for D = 5, all three of the three available representations 
contain a singlet. 
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FIG. 1: Graph of Mod(su(3), n) vs. n, with only points n = 1 (mod 3) shown, fitted by the curve 
{a/n) exp{bn^) with parameter values given by Eqn. 0. The points with n = (mod 3) could 
also be fitted to a smooth curve, as could those congruent to 2 (mod 3); those curves have slightly 
different parameter values. 
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FIG. 2: Fractions of su(3) representations that contain a singlet vs. D, for D < 110. There are 
three distinct curves corresponding to the distinct congruence classes of D modulo 3. 
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Define Modi(su3, D) by the same formula as (0) but witli tlie sum restricted to partitions 
containing 1 as one of tfie elements; tlie fraction of SU3 representations that contain a singlet 
is then Modi/Mod. In Fig. |21 we plot Modi/Mod as a function of D. The creation of 
this plot could not have been possible without the computational speedups suggested by 
Theorems Q and |2l which is most likely why it has never appeared in the literature before. 

Fig. 121 has many fascinating features, which should be of interest to pure mathematicians 
as well as to physicists. Beginning around D = 10, one can notice three series in the 
plot which seem to converge. These three series correspond respectively to the cases D = 
0,1,2 (mod 3). 

Within each series, the points are very regular and seem to line themselves up along a 
smooth curve, which at first seems mysterious. It's possible that an explanation for this 
behavior may be provided by some generalization of Conjecture ^ If we fix the residue class 
of n modulo 3, and then find that 

Modi(3U3,ra) ~ (cLi/n) exp{g{n)) 

for some function g, then we have 

Modi / Mod ~ — exp {bn" - g{n)) . 

If g{n) is such that hn'^ — gin) is negative but approaches zero monotonically from below, 
then we obtain the behavior observed in Fig. |21 



B. Noiseless Subsystems 

Decoherence-free subspaces do not provide the most general method for decoherence-free 
encoding of quantum information. Knill, Laflamme, and Viola discovered a method 
for decoherence-free encoding using subsystems. Zanardi soon thereafter realized that this 
allowed a unification of many seemingly unrelated ideas for reducing decoherence jisl ]. 
Kempe et al. developed a general theory of universal quantum computation based on the 
decoherence-free (or noiseless) subsystem concept |22|. 

Irreducible representations of the error algebra do not possess noiseless subsystems. We 
therefore devote this section to determining how many representations contain noiseless 
subsystems of various sizes. Our results are presented in Figs. El and El and the surrounding 
discussion. 
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A representation of g lifts to a unique associative algebra homomorphism (/> of the 
universal enveloping algebra U{q), by the universal property most elegantly expressed in the 
commutative diagram 

0-^W(g) (8) 




End(V") 

The action of (f) is simply to convert the tensor product to matrix multiplication, i.e. (j){x (S> y) 
= (j){x) ■ (f){y), etc. 

Lemma 2. If(f) is an irreducible faithful representation and if g is a semisimple Lie algebra, 
then (f) is surjective. 

(The proof of lemma El is an easy exercise.) 

Each irreducible representation (p will have the property that (f){U{g)) will be the full 
matrix algebra End(V). On the other hand, noiseless subsystems take advantage of a basis 
in which (f){U{g)) is block diagonal. Therefore, these blocks correspond to a decomposition of 
the Hilbert space Ti into irreducible subrepresentations of the original representation coming 
from the F-operators, 

ni rir 
H=(0K.,)©---©(0F,.) 

where the subspaces are labelled by dimension, dim(VrfJ = di. In this situation, we can send 
an iV-dimensional vector, where = Yli'^^i' through the quantum channel unaffected by 
decoherence. 

If the representation of the F-operators does indeed have such a decomposition, then it 
corresponds naturally to a partition 

r 

D = ^ Uidi 
1=1 

One may now compute the percent of all representations in dimension D which preserve an 
A^- dimensional vector, for various N < D. Recall that R^, which was defined in denotes 
the set of positive integers which are the dimension of an irreducible SU3 module. 

Theorem 3. Let each partition p of D be expressed as D = Yl'i=i^iPi- number of 

representations in dimension D which preserve an N -dimensional vector is given by 

^rit + i{p^) - 1' 



E n 

PG-P(i?3,-D) *=1 
with N=^^ rii 



n.: 
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FIG. 3: the fraction fdiN) of all representations in dimension d = 76 which preserve an A^- 
dimensional noiseless subsystem, as a function of N. The peak is at ~ d/5. 

where the sum is over partitions p with N total repetitions, and with each pi & R^. 

Like theorem^ this theorem also gives a method for computing something which would 
otherwise have seemed elusive. In this case we can plot, for any particular dimension d, 
the fraction fd{N) of all representations in dimension d which preserve an A^-dimensional 
noiseless subsystem, as a function of A^. The results (figs. El and IH) are illuminating. 

First note from the definition of fd{n) that 

d 

n=l 

since every representation is counted once. Therefore, we may view fd{n) as the probability 
distribution function for a random variable. Further, if c? 3> 1, then the collection of points 
{(ra, fd{n)) : n = 1 . . .d} forms a smooth curve, ^ which must then subtend unit total area; 
see figure El 

Figure El shows the exact values of /lool'^) together with a fit to the standard inverted 
beta statistical distribution, defined by 



(9) 



^ For SU3, it is sufficient to take d > 50 to see the smooth shape. 
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FIG. 4: Data points represent exact computations of /ioo(-^)) the fraction of representations in 
dimension d = 100 which preserve an A^-dimensional noiseless subsystem. The points are fitted to 
an inverted beta distribution (see equation @) with a ~ 116.907 and [5 ~ 6.091. Comparing with 
figure 01 which has d = 76, we note that in both cases, the peak is near d/b. 

where a, /3 > and ria^p = x"^^{l — x)^^^dx is a normalization constant. 

Since the distribution fd{n) is not mathematically known to take the form 0, it is of 
interest to determine how close the agreement is. It is important that the points in Figure 
0] are not data; they are exact values for fractions. Therefore, the relevant "scale" in the 
plot is the dimensionless number 1, and a particularly meaningful measure of accuracy is 
simply Af := (l/d) Yl'i=i IVi ~ the average deviation of the model from the "data," 

measured with respect to the only meaningful "scale," which is 1. For the example depicted 
in Figure m we compute = 5.2 x 10~^. 

This agreement is remarkable, and seems to improve as d increases; therefore, we come to 
the unexpected conclusion that fd{n) is asymptotically (for d ^ 1) equivalent to an inverted 
beta distribution! 

Further, inspection of many examples shows that for sufficiently large d, the shape of the 
function fd{n) does not depend strongly on d. More precisely, as d is increased, the function 
receives overall scale factors for the horizontal and vertical axes. In all observed examples, 
the curve always has a peak at around d/5. Given that the distribution function fd{n) is 

13 




completely insensitive to the details of the system Hamiltonian Hs and to the details of 
the heat bath to which the system is coupled, and displays a simple scaling behavior with 
respect to the dimension of the system Hilbert space, one is tempted to term this "universal 
behavior." 

IV. CONCLUSIONS AND OUTLOOK 

As soon has one has appreciated the importance of representations of the sUn Lie algebras 
in physics, natural questions arise. How many of these representations exist, up to equiva- 
lence, for fixed d = dim(7i)? How does the number grow with rf? Suppose that we reduce 
each representation into irreducibles. What is the statistical distribution function which 
describes the relative frequencies of the different possible reductions? The present work has 
given answers (in some cases only partial) to each of these questions, while interpreting the 
results within the unified understanding of quantum decoherence that has emerged in the 
last decade [3, 13- Many unanswered questions remain, such as a proof of Conjecture HJ 
a better explanation of Fig. |21 than the heuristic one given in Sec. IIII Al analogues of these 
results for other Lie algebras, and a fundamental derivation of the "(i/5 rule," by which 
we mean the curious fact that the curves plotted in Sec. ITTTbI seem to always peak around 
N = d/5. 
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APPENDIX A: COMPUTING THE NUMBER OF IRREDUCIBLE REPRESEN- 
TATIONS IN ANY DIMENSION 

In this appendix, we describe how to efficiently compute the function ^{n) that was used 
elsewhere in the paper. 

For a Young diagram A, we let d{X,g) denote the dimension of the corresponding irre- 
ducible module of the Lie algebra q. Unless noted otherwise, we will henceforth assume 
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Q = sUr, and denote the dimension by d{X, r). Further, we often set r = 3, in which case we 
will suppress r from the notation. Young diagrams for SU3 are characterized by row lengths 
Hi, n2 with ni > and < ^2 < rii. The fundamental representation is A = (1, 0) while the 
adjoint is A = (2,1). The representations {ni,n2) and {ni,ni — 722) are conjugate to each 
other and have the same dimension. Representations of the form (2n, n) are self-conjugate. 



The dimension given by the Weyl character formula 



2J is then 



d{ni,n2) = ^{rii + 2){n2 + l)(rai - ^2 + 1) 
= xy{x - y)/2 

where x = ni + 2, and y = n2 + 1. 

We now compute the total number Nij of irreducible representations of 5U3 with dimension 
less than D, by finding, for each fixed rii, the number of n2 which satisfy d{ni,n2) < D, 
and then sum over rii. Expressing the sum in terms of x = rii + 2 gives simpler notation. 
We give the result as theorem EJ omitting the lengthy but straightforward proof. For a real 
number 7 G M, we let [7J denote the greatest integer less than 7. Similarly, [7] denotes the 
least integer above 7. 

Theorem 4. Let No denote the total number of irreducible representations of su(3) with 
dimension less than D. Then No is given exactly by 

Id 

^jr^f/n — — 



Nd = Iid{Id-1)- ^ (M -[!/-]+ !)• (Al) 



where ko := [2\^\, Id := L|(l + VTT8D)\, and 

y^ = ^-^±l{x'-8D/xf\ 

Corollary 1. The exact number ^{D) of irreducible representations in dimension D is given 
by 

i{D) = Nj,^.-Nj,.. (A2) 



For determining the number of irreducible representations of dimension eqns. ()A1|) - 
21) provide a radical computational speedup over the naive algorithm of enumerating all 
possible Young diagrams and computing the dimension of each. These equations may be 
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implemented with an optimized C program. 
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